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ABSTRACT 

The Carnegie-Irvine Galaxy Survey (CGS) is a comprehensive investigation of the physical properties of a com- 
plete, representative sample of 605 bright {Bt < 12.9 mag) galaxies in the southern hemisphere. This contribution 
describes the isophotal analysis of the broadband (BVRI) optical imaging component of the project. We pay close 
attention to sky subtraction, which is particularly challenging for some of the large galaxies in our sample. Ex- 
tensive crosschecks with internal and external data confirm that our calibration and sky subtraction techniques are 
robust with respect to the quoted measurement uncertainties. We present a uniform catalog of one-dimensional 
radial profiles of surface brightness and geometric parameters, as well as integrated colors and color gradients. 
Composite profiles highlight the tremendous diversity of brightness distributions found in disk galaxies and their 
dependence on Hubble type. A significant fraction of SO and spiral galaxies exhibit non-exponential profiles in 
their outer regions. We perform Fourier decomposition of the isophotes to quantify non-axisymmetric deviations 
in the light distribution. We use the geometric parameters, in conjunction with the amplitude and phase of the 
m = 2 Fourier mode, to identify bars and quantify their size and strength. Spiral arm strengths are characterized 
using the m = 2 Fourier profiles and structure maps. Finally, we utilize the information encoded in the m = 1 
Fourier profiles to measure disk lopsidedness. The databases assembled here and in Paper I lay the foundation for 
forthcoming scientific applications of CGS. 

Subject headings: atlases — galaxies: fundamental parameters — galaxies: general — galaxies: photometry — 
galaxies: structure — surveys 



1. INTRODUCTION 

This paper, the second in a series, presents the isophotal anal- 
ysis for the optical images of the Carnegie-Irvine Galaxy Sur- 
vey (CGS), a detailed study of a statistically complete sample 
of nearby, bright galaxies in the southern sky (Ho et al. 2011, 
hereinafter Paper I). The immediate aim of this paper is to re- 
duce our extensive set of images to a uniform database of one- 
dimensional (1-D) radial profiles of surface brightness and geo- 
metric parameters, on which much of our subsequent scientific 
analysis will depend. Although we intend to apply more so- 
phisticated methods of analysis to the images (Peng et al. 2010; 
S. Huang et al. 2011, in preparation), the 1-D analysis already 
contains a wealth of useful information that can be exploited for 
science. Moreover, 1-D analysis has the virtue of simplicity. It 
can be efficiently applied to a large sample of objects, allowing 
a quick overview of the global properties of the survey. 

The brightness profiles of galaxies have long helped to guide 
our understanding of their physical nature. Despite the visual 
complexity of their images, the 1-D radial brightness profiles 
of galaxies in the nearby universe actually show a surprising 
degree of order. De Vaucouleurs (1948) first noticed that the 
light distributions of elliptical galaxies generally follow a r'/"* 
profile, which has been inter preted as a signatu r e of dissipa - 



since been generally adopted as the standa rd formula for fitting 
the brightness profiles of ellipticals (e.g., lOraham etal] 1 19961 : 
iTruiillo et aLll200ll: iKormendv et alj|2009[) . Our modern view 
of bulges has also grown steadily more complex over time. 
Once thought to be mini-ellipticals with r'/"* profiles, bulges, 
too, are now known to be better described by a Sersic r'^" 
function (lAndredakis & Sanders 199 4 ; Andredakis et al. 19951 : 
Ide JongI [19961: ICourteau et al.1 119961: fMacArthur et al.l l2003b . 



The Sersic indices of bulges have a br oad distribution of ob- 
served values, from n < 1 to » > 4 (e .g., [Mac Arthur et al.ll2003l : 
iF isher & Drorv| l2008': Gadottil |2008h . and it is argued that they 
reflect different formatio n physics. Spheroids w ith « ^ 2 are re- 
garded as pseudobulges jFisher & Drorvl[2008l) . which formed 
through internal, secular processes, while those with « ^ 2 
are classical bulges, which, like the ellipticals, were assem- 
bled more rapidly, most likely with the assistance of mergers 
(IKormendv & Kennicut3l20o4 . 

The brightness profiles of the disks of SO and spiral galax- 
ies have been traditionally described by a single exponen- 
tial function (de Vaucouleurs 1959; Freeman 1970), which 
arises as a natural consequence of v iscous transport in a 
disk (lYoshii&Sommer-Larsenl 119891: IZhang & Wvsd 120001: 



tionless formation processes jvan Albadalll982t iKatd Il991h. 



Later studies, beginning with ICaon et ak r (ll993h . increasingly 
recognized that many ellipticals, in fact, do not strictly follow 
the r'/'* law, but inst ead are better described by the more gen- 
eral r'/" function of lSersid ( Il968h . of which de Vaucouleurs' 
law is a special case (« = 4). Indeed, the Sersic function has 



[Ferguson & Clarke[ 2001 : Slvz et alj2002 ). p erhaps mediated 
by st ar format ion and feed back processes dRobertson et al.l 
[20041 : iGover nato et al.ll2007[) . In actuality, very few disks are 
so simple. Many possess breaks and inflections in their outer 
radial profile dvan der KruitI 19791: van der Kruit & Searlell981l: 
Matth ews & Gana gheiill997l:lPohlen et al.ll2000l; Ide Grijs et al.l 
2OOTI: lErwinll2005UErwin et al.ll2005l: iPohlen & Truiilloll2006l) . 



No general consensus yet exists as to their cause, but they offer 
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Fig. 1. — Illustration of how we determine the sky radius for the i?-band image of NGC 1400. Left: star-cleaned image, showing the full field-of-view of 8f9x8f 9. 
Right: radial profile of the isophotal intensity, in units of ADU pixel"' . The veilical dash-dotted line marks the radius where the isophotal intensities start to oscillate 
rather than continue to decrease; this radius, designated by the keyword SKY_RAD in the image header, denotes the region outside of which the sky dominates. The 
con'esponding isophotal ellipse is overplotted in the left-hand panel. 



important clues t o a host of physical processes pertinent to 
galaxy formation ( I van der Kruit & Freemanll201 lb . 

Apart from intensity profiles, isophotal analysis of galaxy 
images yields other useful diagnostics. The radial varia- 
tion of the ellipticity and position angle, for example, pro- 
vides an efficient means to identify bars and to quantify 
their length and strength (e g., Laine et all 120021: lErwinI [20051: 
iMenendez-Delmestre et al.l |2007|) . Fourier decomposition of 
the isophotes provides yet another method to probe non- 
axisymmetric perturbations in the light distribution. The rel- 
ative amplitude of the m = 2 mode, in combination with 
its phas e angle, has been shown to be effective in isolat- 
ing ba rs dElm egreen & Elmegreen 1985VButa"1986';'Ohta et al.l 
19901) and spira ls (Elmegreen et al. 1989; Rix & Zaritsky 19951 



Odewahn et alj l2002h . Both of these features are common 



ble integrated colors and color gradients (Section 6). The prod- 
ucts from the isophotal and Fourier analysis are used to quan- 
tify the strengths of bars (Section 7), spiral arms (Section 8), 
and lopsidedness (Section 9). Section 10 assesses the reliabil- 
ity of our measurements using internal and external tests. Sec- 
tion 1 1 gives a brief summary and an outline of future plans. 
The database of isophotal parameters is described in the Ap- 
pendix. 

2. SAMPLE PROPERTIES 

The CGS covers a statistically complete sample of 605 
bright, nearby galaxies of all morphological types in the south- 
ern hemisphere, with B-band total magnitude Bj < 12.9 and 
5 <Q°. These very general selection criteria enable us to probe 
galaxies with a broad range of physical properties and mor- 
phologies. The primary parent sampl^ comprises 17% ellip- 
ticals, 18% SO and SO/a, 64% spirals, and 1% irregulars. The 
bulk of the sample is relatively nearby (median Dl = 24.9 Mpc), 
luminous (median Mg^ = -20.2 mag), and well resolved. The 
typical seeing of CGS is ~ 1", and the sample has a median 
isophotal angular diameter of D25 = 3. '3 at a surface brightness 
level of = 25 mag arcsec"^. 

Paper I describes the observing strategy, data reductions, and 
photometric calibration of the optical imaging component of 
the project. We only repeat a few essential details here. The 
broadband BVRI images have a field-of-view of 8.'9x8f9 and 
a pixel scale of 0"259, which is well matched to the good see- 
ing typically achieved with the du Pont 2.5-m telescope at Las 
Campanas Observatory. The median seeing of the survey, as de- 
termined from over 6000 science images, is 1 f'17, 1 "1 1, 1 "01, 
and 0"96 in the B, V, R, and / band, respectively. A little more 
than half of the galaxies were observed under photometric con- 
ditions, with median photometric errors of 0.08, 0.04, 0.03, and 

^As described in Paper 1, we observed an additional 1 1 galaxies that do not formally meet the selection criteiia of CGS. We still analyze them here but will not use 
their results to draw statistical inferences on the sample. 



constituents in disk galaxies, and both are thought to play 
a dynamical role in facilitating angular momentum transport 
and driving secular evolution. Likewise, a significant frac- 
tion of disk galaxi es exhibits global lops i dedness in their stellar 
light distribution jRix & ZaritskvllT995t IZaritskv & Rixlll997t 
iBournaud et al.l 12005 : Reich ard et al.ll2008l), whose main cul- 
prit remains in disput eHog & Combes 2008). As shown by 
iRix & Zaritskvl (11995!), this type of non-axisymmetric pertur- 
bation is again conveniently revealed through Fourier analysis 
of the isophotes, in this case through the m=l mode. 

This paper is organized as follows. Section 2 gives a brief 
overview of the CGS sample, the observations, and some basic 
characteristics of the images. Section 3 describes our method 
of sky subtraction. The procedural details of isophotal analysis 
are presented in Section 4, including our method for extracting 
geometric parameters, surface brightness profiles, and Fourier 
components. We generate composite light distributions (Sec- 
tion 5) to identify statistical trends in disk profiles, and assem- 



Fig. 2.— Left: B-band image of NGC 1400, binned by 20 X 20, with the ellipse in Figure [T]overplotted. We show the full field-of-view of 8;9x8f9. The 
sky value and its uncertainty are simply the mean and standard deviation of the pixel values outside of this ellipse, after excluding the masked objects, which are 
shown as black regions. Right: normalized histograms of the background pixels of the original and the binned images; their peak positions are 0.3060 ± 0.0102 and 
0.3059 ±0.0011 ADUs"' pixel" ', respectively. (A color version of this figure is available in the onhne journal.) 



0.04 mag for the B, V, R and / filters, respectively. We devised 
a calibration strategy to establish an approximate photometric 
zero point for the non-photometric observations, for which the 
corresponding photometric errors are 0.21, 0.1 1, 0.10, and 0.09 
mag. After correcting for large-scale gradients in the back- 
ground, the flatness of the final images is about 0.6%, and the 
typical depth of the surface brightness, defined as 1 cr above the 
background, has a median value of yu « 27.5,26.9,26.4, and 
25.3 mag arcsec"- in the B, V, R, and / bands, respectively. 

We derived a number of data products from the reduced, 
calibrated images. These include red-green-blue color com- 
posites generated from the B, V, and / bands, images cleaned 
of foreground stars and background galaxies, a stacked image 
from a weighted combination of the four filters optimized to 
enhance regions of low surface brightness, structure maps de- 
signed to accentuate high-spatial frequency features, and color 
index maps from different combinations of the filters. 

3. SKY DETERMINATION 

Sky determination is a crucial, fundamental step in the data 
analysis. Many of the basic galaxy parameters we are inter- 
ested in measuring (magnitudes, colors, characteristic size and 
brightness level, etc.) are predicated on having the sky level 
properly subtracted. Importantly, under-subtraction or over- 
subtraction of the sky value can introduce spurious curvature 
into the brightness profile, especially i n the faint, outer regions 
of the galaxy (e.g.. lErwin et al...2008l) . iMacArfliur et alj (l2003h 
studied the influence of the sky value on the bulge and disk pa- 
rameters for a sample of spirals, and concluded that the disk, 
but to a lesser extent even the bulge, parameters are sensitive to 
the sky value. 

There are a variety of ways to measure the sky value 
of a CCD image. Science data such as ours, however, 
wherein an extended object fills a substantial portion of the 
chip, pose unique challenges. This is especially so be- 



cause the background of our images is not always entirely 
uniform (Paper I). We adop t a tw o-step approach. As in 
iNoordermeer & van der Hulsa ( |2007|) . we generate the isopho- 
tal intensity radial profile of the galaxy to the edge of the field 
to determine the radius beyond which the sky background dom- 
inates the signal. As illustrated in Figure [T] the transition from 
the galaxy's outer boundary to the sky-dominated region man- 
ifests itself as flattening of the radial profile, beyond which 
it oscillates about a constant intensity level. To be specific, 
we define the outer radius to be the first data point where the 
measured isophotal intensity rises instead of decreases mono- 
tonically. Typically the outer radius is large enough to avoid 
the spiral arms or other features that may cause a real rise 
in the outer brightness profile. In t he standard procedure of 
INoordermeer & van der Hulsfl ( |2007|) . the average value of this 
isophotal intensity outside the outer radius and the associ- 
ated standard deviation gives estimates of the sky level and 
its uncertainty. However, this technique is relia ble only if th e 
background is uniform and well measur ed (Erwin et alJliOOSh . 
The field-of-view of our images is typically only a factor of 
^2 larger than the galaxies, generally too marginal to pro- 
vide enough data points in the sky-dominated region to yield 
a statistically robust measurement of the background and its 
error The situation is further exacerbated by the occasional 
presence of residual large-scale non-uniformities in the back- 
ground. In view of these complications, we use Noordermeer 
& van der Hulst's method only to determine the radius of the 
sky-dominated region (Figure[T]i, which we record in the image 
header under the keyword SKY_RAD. 

To estimate the actual sky value and its associated erro r, we 
follow a method similar to that used bv lErwin et al. I (l2008h . We 
first smooth the original 2042 x 2042 pixel image by binning it 
down to a 102 x 102 pixel image. This highlights underlying 
large-scale, systematic fluctuations in the background, which is 
the main factor that ultimately limits the accuracy with which 
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Fig. 3. — Empirical relationship between skygi/skyreai and Rat/Ris, for 
each of the four filters. For Rfn/Ris > 1-3, skyj), fa sky^jj^ji as indicated by the 
best-fitted red solid lines with slopes close to 0, but when Rfn/R25 < 1.3, sky^, 
is systematically larger than sky^^j,!. The green solid lines represent the best-fit 
linear relations given by Equations (l)-(4). (A color version of this figure is 
available in the online journal.) 

we can determine the sky level, and hence the final sensitiv- 
ity of the surface brightness of our images. The value of each 
binned pixel is the mean value of all the pixels inside a 20 x 20 
pixel box, after excluding field stars and background galaxies 
identified in the image's object mask (Paper I). The background 
pixels, then, are defined to be all the pixels in the binned image 
outside of the isophotal ellipse marked by SKY_RAD. The sky 
level is simply the mean of the background pixel values, and 
its uncertainty is the standard deviation of the pixel flux distri- 
bution; these values are stored in the image header under the 
keywords SKY_VAL and SKY_ERR. (The standard deviation 
of the mean, useful for other applications, is stored separately 
under the keyword SKY_SIG.) We have opted to compute the 
mean of the pixel flux distribution, rather than its median or 
mode, but in practice this makes little difference (they gener- 
ally agree to ^0.1%) because the shape of the distribution is 
highly symmetrical. This reflects the robustness of our estimate 
of the sky-dominated region and the effectiveness of our ob- 
ject masks in rejecting faint halos around foreground stars and 
background galaxies. 

We have tested the effect of choosing different scales for the 
smoothing, varying the binning box sizes from 5 to 50 pixels. 
While the average sky value remains stable, the width of the 
pixel distribution decreases with increasing smoothing length, 
leveling off to a near-constant value for box sizes ^20 pixels 
(~5")- We interpret this to represent the typical scale of large- 
scale systematic fluctuations in the sky background. We choose 
a box size of 20 x 20 pixel as a reasonable compromise in or- 
der to retain a statistically significant number of data points to 
compute their average and standard deviation. 

Figure[2]illustrates our method of sky estimation, using a B- 
band image of NGC 1400. The image has been binned 20 x 20, 
and the foreground stars and background galaxies have been 
masked out. The sky value is the mean of pixel values outside of 



the ellipse, after excluding the masked objects, and the error is 
simply the standard deviation of the sky pixels in the binned im- 
age. The right-hand panel shows normalized histograms of the 
sky pixels of the original and binned image. Clearly they peak 
at nearly identical locations (the peaks of the two histograms 
differ by 0.0001 ADU s"' pixel"'), but the distribution for the 
binned image is narrower than that of the original image by a 
factor of 10. 

The above-described strategy for sky determination can only 
be applied to galaxies with angular diameters D25 <; 5'-6'. For 
the ^15% of the sample more extended than this, it is difficult 
to impossible to determine the radius of the sky-dominated re- 
gion and obtain robust statistics for the sky pixels, and we must 
resort to a more indirect approach. The signal in the outer re- 
gions of the CCD frame consists of galaxy light plus a constant 
sky background. Assuming that the galaxy component can be 
modeled by a single Sersic function, we can fit the observed 
light profile to solve for the underlying sky value. We perform 
the fitting on the 1-D surface brightness profile (Section 4), af- 
ter excluding the central regions of the galaxy and other features 
such as the bulge, the bar, or strong spiral arms, if present. Sim- 
ple experimentation shows that the best-fit sky value depends 
on the fitting radius relative to the size of the galaxy. Clearly, 
if the fitting radius is large compared to the outer edge of the 
galaxy, the sky will be well determined; however, if the fitting 
radius lies substantially interior to the main body of the galaxy, 
the inferred sky value will depend critically on how well the 
Sersic model represents the intrinsic light profile of the galaxy. 

We devise an empirical correction, as follows. We select sev- 
eral galaxies that (1) have relatively simple structures, (2) are 
small compared to the CCD's field-of-view, and (3) have well- 
determined sky values. Then, we fit their surface brightness 
profiles with different fitting radii (Rfn), to mimic the actual 
situation in galaxies that are too angularly extended to have 
a reliable sky determination. The resulting fitted sky value, 
skyjj,, is then compared with the independently known, correct 
value sky Figure |3] shows sky^j/sky^g^, vs. /?fit//?25, where 
R25 = 0.5Z)25. We can see that skyfjj/skyj.g^] ?» 1 when Rfn/Ris 
^ 1.3. When Rfn/Ris < 1.3, skyjj^ overestimates sky^^^,, but 
it does so systematically, in such a way that we can apply an 
approximate empirical correction to recover the true sky value. 
The best-fitting linear relations (and their associated rms scat- 
ter), are s follows. 



Bband 

sky fit 



Sky,-eal 

Vband 

Sky^eal 

7? band 

skyfit . 

Sky^eal 

/band 

skyfi, 



= 1.364-0.327 X 



^fit 



(7 = 0.05, 



(1) 



1.427 



1.309 



-0.389 x:^, 

^25 



-0.287 x:^, 

^25 



(7 = 0.033, (2) 



(7 = 0.042, (3) 



sky, 



= 1.070-0.060 X 



real 



R^5' 



(7 = 0.020. 



(4) 



The best-fit sky value, skyj^j, and its associated statistical er- 
ror, (Tfit, are recorded under the header keywords SKY_VAL and 
SKY_ERR with the comment "fitted sky value." If the above 
empirical correction to sky^j is necessary, we fold the scatter of 
the correction relation into the error budget. 
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4. ISOPHOTAL ANALYSIS 

4. 1 . Geometric Parameters and Surface Brightness Profile 

The IRAF0 task ellipse is commonly used to measure the 
surface brightness profiles of galaxies ( e 
19941; iMilyang- Jensen & J0rgensenl 119991; 



Jogee et al.ll2004lAguerri et al.ll2005 



Noordermeer & van der Hu lst 2007; 



Silva & ElstonI 
Laine e t al 1 120021; 



Marinova & J ogee 200J7 
Barazza et al. 2008) 



Fol- 



lowing the iterative method of Jedrzejewski (1987), we fit the 
isophotes of the galaxy with a set of ellipses. This is motivated 
by the fact that the isophotes of most galaxies, especially early- 
type systems such as ellipticals and lenticulars, are quite close 
to ellipse^ In our implementation, the ellipses are sampled 
along the semi-major axis of the galaxy in logarithmic inter- 
vals, starting from r « 0"3 and increasing the radius of each 
successive ellipse by a factor of 1.1. After reaching the out- 
ermost ellipse, the fitting reverses direction and moves toward 
the galaxy center from r « 0"3, with each subsequent radius 
decreasi ng by a factor of 1 . 1 . 

As in iNoordermeer & van der Hu"is3 (|2007|) . we determine 
the isophotal geometric parameters of the galaxy in two steps. 
In the first step, we estimate the center of the galaxy. Ellipses 
are fitted to the /-band image with the center, position angle 
(PA), and ellipticity (e) set as free parameters. We use the I- 
band image as the fiducial reference because of its relative in- 
sensitivity to dust extinction and young stars, and because it 
generally has the best seeing. The center of the galaxy is the 
average central position of the ellipses inside ~5"-7". In im- 
ages of regular galaxies the center of the best-fit ellipses often 
converge to a well-defined value, with rms « 0."015. However, 
in galaxies with dusty nuclear regions, the best-fitting central 
isophotes may give a poor measure of the true center. In such 
cases, a better estimate of the true center comes from isophotes 
at intermediate radii, ■~10"-30", far enough to be undisturbed 
by central dust but yet sufficiently close to the nucleus to give a 
faithful measure of its position. The typical uncertainty of the 
central positions estimated in this way is '-^0"02. 

Next, we fix the center just determined and run ellipse again, 
while still setting e and PA free. Our goal is to determine 
the characteristic e and PA of the galaxy based on its best-fit 
isophotes. Typically we take the average value of these param- 
eters in the outer regions of the galaxy, where the intensity is 
about I a above the sky, as their characteristic values and use 
their standard deviations over that region as the uncertainties. 
These parameters usually converge to a constant value within 
that region, with variations of ~0.04 for e and ^2° for PA. 
However, the intrinsic geometric parameters of some galaxies 
can be distorted by mergers or interactions, causing e and PA 
to diverge at large radii. In these cases we simply estimate their 
values manually from the visually best-fitting isophotes near the 
edge of the galaxy. The B, V, R, and / images have their center, 
e, and PA values determined independently, and they are stored 
in their corresponding image headers. 

During this second step, we also record the deviations of 
the isophotes from perfect ellipses, which, as described in 
IJedrzeiewskil (1 19871) . are parameterized by the third (A3,B3) 
and fourth (A4,B4) harmonics of the intensity distribution. The 



A3 and £3 parameters g ive "egg-shaped" o r " heart-shaped" 
isopho tes dCarteJ 119781; IJedrzeiewskil Il987l) . iPeletier et al] 
(1 19901) point out that A3 and B3 appear to be sensitive diagnos- 
tics of dust features in elliptical galaxies. The most interesting 
parameter among them is B4: if it is positive, the underlying 
isophote is disky with respect to a perfect ellip se; a negative 
B4 cor responds to a boxy isophote. Figure 1 of IPeletier et al.l 
11990*) gives examples of the different isophotal shapes for dif- 
ferent values of B3,A4, and B4. 

We run ellipse for a third and final time to extract the av- 
erage intensity of the isophotes, fixing the geometri c param- 
eters to the values determined above (e.g ., Pohlen & Truiillol 
I2OO6I; INoordermeer & van der Hulsj|2007l) . For this step, we 
do not allow the geometric parameters to vary in order to re- 
duce the influence of bars and other non-axisymmetric features 
on the average intensity profile, as well as to reach convergence 
in regions where the signal-to-noise ratio (S/N) is marginal 
dErwin et al.ll2008l) . For consistency, we apply this isophote 
measurement to all the galaxies in our sample. For the lopsided 
galaxies, we also experimented with allowing the isophotal cen- 
ters to be left as free parameters. We find that the typical dif- 
ference in the brightness profile, compared with the fits based 
on fixed isophotal centers, is ^ 0.2 mag arcsec"^. Moreover, 
our tests show that the relative amplitudes of the Fourier terms 
(Section 4.2) decrease significantly when the isophotal centers 
are allowed to be free, due to the fact that the free-fitting el- 
lipses can better trace the distorted disk in the outer part of the 
galaxy to produce very small fluctuations along each isophote. 
This will make the Fourier analysis less effective for detecting 
and quantifying the properties of bars or lopsided structures. 

After subtracting the sky background from the image, the 
surface brightness is calculated from 



/X = -2.5 log 



fexp X A 



+ zpt, 



(5) 



where Iiso is the isophotal intensity after subtracting the sky, fgxp 
is the exposure time of the image in units of seconds, A is the 
pixel area in units of arcsec", and zpt is the photometric zero 
point of the image in units of magnitudes. We calculate the sur- 
face brightness only from those isophotes whose intensities are 
larger than /sky + <Jsky We propagate the errors on /iso into errors 
on fi in magnitude units. The surface brightness profiles in the 
B,V, and R bands are constrained to have the same geometric 
parameters as determined in the / band. 

To construct 1-D color profiles, we blur all the images to a 
common seeing. This is done by convolving the image having 
the better seeing with a two-dimensional (2-D) Gaussian func- 
tion whose full width at half maximum (FWHM) is the quadra- 
ture difference between the two seeing values. The measured 
isophotal ellipses of the unblurred /-band image are used to di- 
rectly calculate the isophotal intensity of the blurred images in 
all the filters. Color profiles follow from straightforward differ- 
encing of one band from another 

4.2. Fourier Analysis 



^IRAF is distributed by tlie National Optical Astronomy Observatory, which is operated by the Association of Universities for Research in Astronomy, Inc., under 
cooperative agreement with the National Science Foundation. 

'An important exception is edge-on galaxies, which are not well suited to ellipse fits. As the bulge and disk components have different profiles and ellipticities, 
their relative contributions change with radius and azimuth (i.e. along major or minor axes). A dust lane, if present, also can strongl y affect the averaged isophotal 
intensity. For e dge-on galaxies it is preferable to extract the isophotal intensities along cuts in the major and minor axis directions (e.g., lde Griisll998l ; lFrvet al.ll999l ; 
IWu et alJ2002i) . 
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Fro. 4. — B-band composite profiles of our sample, divided by morphological type. The elliptical galaxies are normalized at a reference radius of r,ef = R20, and 
the profiles are plotted vs. (r/fjef)'^*- The thick black line corresponds to a de Vaucouleurs r''''' law. The profiles for the disk galaxies are scaled according to the 
scale length h of the disk outside rief . Profiles of Type I, II, and III are marked in red, purple, and blue, respectively. The thick black line corresponds to a pure 
exponential function. (A color version of this figure is available in the online journal.) 



We study the harmonic components of the intensity distribu- 
tion of the isophotes. In our work, we decompose the intensity 
distribution along each ellipse into a Fourier series of the form 

00 

m = Io + Y,tjCosji9 + <j>j), (6) 
>=1 

where / is the intensity (in units of ADU s"' pixel"') on the 
ellipse in the direction 9, Iq is the average intensity of the el- 



lipse, Ij measures the strength of the yth mode in the series, 
and (pj is the corresponding phase angle of that mode. The 
angle 6 is defined to be 0° along the positive y-axis and in- 
creases counterclockwise; 0y = 0° along the positive y-axis and 
increases clockwise. A high S/N is required to derive a statisti- 
cally significant measurement of the high-order Fourier terms 
dNoordermeer & van der Hulsil2007h : thus, we only perform 
this decomposition inside the radius where the average intensity 
of the isophote is 3 dsky larger than the determined sky value. 
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2 3 4 




Sc-Sd 





Fig. 5 . — /-band composite profiles of our sample. See Figure|4]for details. (A color version of this figure is available in the onhne journal.) 



The relative amplitude of the m = 1 (/i//o) and m = 2 (h/Io) 
mode will be especially useful in our analysis, as the y reflect 
the lopsidedness of the galaxy (Rix & Zaritskvl ll995l) and the 
st rength of the bar or spiral arms (.Buta.1986.) , respectively. 

iButal (Il986l) performed a Fourier analysis to study the az- 
imuthal variations of the light distribution of NGC 1433, which 
is also included in the CGS sample. We found good agreement 
between Buta's values of the relative amplitudes of the m = 1 
and 2 modes and those calculated by us. This helps to confirm 
the robustness of our method. 



4.3. Database of 1-D Profiles 

The full database of isophotal parameters for the 605 
galaxies in CGS (including the 11 extras not formally 
part of the survey) is given in the Appendix, as Fig- 
ures 19.1-19.616, as well as on the project Web site 
|http : //cgs . obs . carnegiescience . edu[ 

5. COMPOSITE PROFILES 

Composite profiles can help to highlight characteristic statis- 
tical trends, as well as to isolate interesting outliers, in a class of 
objects. We normahze the surface brightness profiles to a com- 
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TABLE 1: Frequency of Profile Types for the Disk Galaxies 



Hubble Type 


T 


Filter 


Number 


Type I 


Type 11 


Type III 


Other 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


SD-Sn /a 


— 3 < T < 0.5 


B 


75 


16.0% 


26.7% 


57.3% 


0.0% 






y 


78 


20.5% 


26.9% 


52.6% 


0.0% 






R 


75 


24.0% 


24.0% 


52.0% 


0.0% 






I 


81 


30.9% 


19.8% 


49.3% 


0.0% 


Sa-Sab 


f) 5 < T < 2 5 


B 


51 


17.6% 


39.2% 


39.2% 


4.0% 






V 


52 


15.4% 


40.4% 


40.4% 


3.8% 






R 


53 


18.9% 


37.7% 


41.5% 


1.9% 






I 


52 


30.8% 


34.6% 


34.6% 


0.0% 


Sb-Sbc 


2.5 < T < 4.5 


B 


132 


12.1% 


63.6% 


24.3% 


0.0% 






V 


134 


20.1% 


53.7% 


24.6% 


1.6% 






R 


135 


22.2% 


49.6% 


26.7% 


1.5% 






I 


135 


27.4% 


45.2% 


27.4% 


0.0% 


Sc-Sd 


4.5 < T < 7.5 


B 


137 


7.3% 


72.3% 


19.7% 


0.7% 






V 


138 


10.1% 


71.0% 


18.1% 


0.8% 






R 


137 


13.1% 


67.2% 


18.2% 


1.5% 






I 


137 


27.7% 


56.2% 


14.6% 


1.5% 


Sdm-Sm 


7.5 < T < 9.5 


B 


17 


23.5% 


35.3% 


41.2% 


0.0% 






V 


17 


23.5% 


35.3% 


41.2% 


0.0% 






R 


16 


31.3% 


37.5% 


31.2% 


0.0% 






I 


17 


29.4% 


29.4% 


41.2% 


0.0% 



NOTE. — Column 1: Hubble type. Column 2: morphological type index. Column 3: filter. 
Column 4: number of objects. Columns 5-7: fraction with profiles of Type 1, 11, and III, respectively. 
Column 8: other complicated profiles. 



mon reference radius. For the elliptical galaxies in our sample, 
we set Tief to /?20, the radius wherein 20% of the total flux is en- 
closed, and plot the profiles as a function of (r/rref)'/'*. In this 
reference frame, a classical de Vaucouleurs r'/"* profile traces a 
straight line. Figure|4](top left panel) illustrates the now-known 
fact that not all ellipticals obey the r'/"* law, but rather are better 
described by a more general Sersic function with r'/". 

For the disk (SO and spirals) galaxies, we set rret to be roughly 
the boundary between the bulge and disk, and we plot the 
brightness profiles as a function of (r-r,ef)//!, where h is the 
scale length of the disk, as determined by fitting an exponential 
function to the profile outside of r^^f. This choice of coordinates 
helps to reveal possible deviations of the disks from a canonical 
exponential profile (n =1), which appears as a straight line. It is 
apparent that the light distributions of the disk very rarely fol- 
low a pure exponential function (Figure |4]i, especially in their 
outer regions. Most show a downward turn compared to a single 
exponential, but not an insignificant number s how an upward 



turn. This phenomenon is well-known (e.g., Phillipps et al.l 
[19911 lErwinll2005l:lErwin etani2005l: iPohlen & TruiiUdl2006h . 



We code the three profile types with different colors, with 
Type I profiles (no break) in red. Type II profiles (downward 
break) in purple, and Type III profiles (upturn) in blue. The 
scatter among the normalized profiles is smaller in the redder 
bands, indicating that dust extinction and young stars have a 
greater effect on the profile shapes in the blue. Figure [5] shows 
the equivalent montage for the / band. 



Table 1 tabulates the frequency of the different profiles for 
each bin of morphological type for the disk galaxies. Oc- 
casionally, there are some galaxies with complicated surface 
brightness profiles, which cannot be classified as any of the 
three standard types listed above. The profiles for such ob- 
jects are listed as "Other" in Table 10. Type II and Type III 
profiles are common in our sample, and their fractions depend 
on the galaxy morphology. Type II profiles occur more fre- 
quently in late-type disk galaxies (~ 70% among Sc-Sd spi- 
rals), whereas Type III profiles are preferentially found in more 
bulge-dominated, earlier-type systems, especially among the 
SO-SO/a class (- 50% -60%). The fraction of galaxies with 
Type I profiles, on the other hand, seems to be roughly con- 
stant, at ^ 20%, across all morphological types; the fraction 
increases systematically toward redder bandpasses, except for 
the Sdm-Sm galaxies, where the fraction seems to be roughly 
constant, although the number of objects is small. A detailed 
analysis of the different profile types and their dependence on 
other physical parameters will be presented in a separate paper. 

6. COLOR INFORMATION 

Table 2 presents integrated colors and color gradients for the 
sample, c orrected for foregroun d Galactic extinction using val- 
ues from lSchlegel et all (Il998h . We Hst B-I,V -I, wA R-I, 
from which other color combinations can be readily derived; 
we use total magnitudes within the last reliable isophote (1 cr 
above the sky), as given in Table 4 of Paper I. For each of 



Since the brightness profile in the outer regions of the galaxy depends sensitively on the accuracy of the sky subtraction, we exclude objects with unreliable sky 
determination. We also omit galaxies whose light distribution is severely adversely affected by very bright foreground stars, by excessively crowded field stars, or by 
an interacting neighbor The excluded objects are flagged in Table 2. Note that a star that is bright and excluded in one filter may not be equally bright or rejected in 
another; therefore, the number of objects in a morphological bin is not the same among all the filters. We further omit the 1 1 extra galaxies that do not formally meet 
the CGS selection criteria. 



TABLE 2: Integrated Colors and Color Gradients for the CGS Sample 



Index 


Name 


B - I 


V -I 


R-I 


V(B - /)in 


V(B - J)out 


V(l^- J)in 


V{V - J)out 


V(i?-/)i„ 


V{R - J)„ut 


Notes 








(mag) 


(mag) 


(mag) 
















(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


1 


ESO 


009-GOlO 


1.83±0.06 


1.22±0.05 


0.68±0.05 


-0.19±0.05 


0.30±0.77 


-0.08±0.04 


0.01±0.72 


0.00±0.04 


0.11±0.73 




2 


ESO 


027-GOOl 


1.31±0.30 


0.90±0.21 


0.52±0.19 


-0.67±0.17 


-0.19±0.73 


-0.32±0.10 


-0.07±0.72 


-0.11±0.08 


0.03±0.73 




3 


ESO 


027-G008 


1.96±0.30 


1.26±0.21 


0.63±0.19 


-0.41±0.08 


0.80±0.66 


-0.18±0.07 


0.69±0.65 


-0.07±0.07 


0.61±0.61 


Extra 


4 


ESO 


056-G115 




















BVRI^ 


5 


ESO 


060-G019 


1.42±0.30 


1.08±0.21 


0.65±0.19 


0.23±0.13 


0.55±1.02 


0.19±0.12 


0.47±1.00 


0.20±0.12 


0.39±0.98 




6 


ESO 


091-G003 


1.67±0.05 


1.04±0.05 


0.49±0.04 


-0.12±0.07 


-0.65±0.36 


0.01±0.06 


-0.25±0.34 


0.08±0.05 


0.00±0.34 


Extra 


7 


ESO 


097-G013 








O.OliO.lO 


-l.lOil.lO 


0.07±0.09 


-0.27±1.09 


0.09±0.08 


-0.39±1.04 


BVRI^' 


8 


ESO 


121-G006 


1.97±0.14 


1.42±0.09 


0.82±0.08 


-0.77±0.14 


0.14±0.55 


-0.40±0.12 


0.38±0.51 


-0.13±0.10 


0.44±0.56 


P 


9 


ESO 


121-G026 


1.68±0.14 


1.15±0.09 


0.63±0.08 


-0.52±0.08 


-0.89±0.29 


-0.19±0.05 


-0.46±0.28 


-0.06±0.04 


-0.19±0.28 




10 


ESO 


136-G012 


1.31±0.39 


1.12±0.26 


0.75±0.19 


-0.30±0.25 


0.06±2.00 


-0.08±0.23 


0.01±1.90 


0.06±0.24 


0.28±2.01 


BVRI^ 



Note. — Column 1: numerical index. Column 2: galaxy name. Columns 3-5: integrated colors derived from the total magnitudes within the last reliable 
isophote (1 a above the sky). Columns 6-11: color gradients of the galaxy interior and exterior to R^q measured in the / band (see Section 6). The 
units are Amag per dex in radius. The colors and color gradients have been corrected for Galactic extinction. Column 12: notes: "extra" = object not 
formally part of the main sample; filters = image of that filter has problem, which is indicated by the table mark. (This table is available in its entirety in 
a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and content.) 

^Galaxy adversely affected by a very bright star. 

''Galaxy adversely affected by very crowded field stars. 

'^Galaxy has unreliable sky value determined from model fitting. 

'^Galaxy distorted by an interacting neighbor. 

''No magnitude zero point available. 
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V{B-I) V{B-I) 

Fig. 6. — Normalized histograms of the inner and outer B — I color gradients, divided by morphological type. The inner and outer color gradients are repre- 
sented in red and blue solid histograms, respectively. The vertical dashed lines in each panel mark the adopted boundaries for negative (V(B — /) < —0.1), flat 
(-0.1 < V(B-/) < 0.1), and positive (V(B-/) > 0.1) color gradients. A positive color gradient means that the color becomes redder outward, while a negative 
value indicates that the color gets bluer outward. (A color version of this figure is available in the online journal.) 
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Fig. 8. — Normalized histograms of tlie inner and outer R-I color gradients. See Figure |6]for details. (A color version of this figure is available in the online 
journal.) 
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these color combinations, we also calculate two simple mea- 
sures of the color gradient, after resampling the color profiles 
with 300 equally spaced data points in linear space, to over- 
come the heavily sampled points in the central region. Sim- 
ilar to iTaylor etaP (l2005h . the color gradient is derived both 
inside and outside the half-light radius, 7?5o, as determined in 
the / band. The inner region ranges from 3 times the seeing 
FWHM (to avoid spurious effects from seeing mismatches) to 
RsQ, while the outer region extends between R^o and 2.5 Rso- 
Tests show that most of the color profiles of the CGS galaxies 
can reach further than 2.5 R50. We confine the color gradient 
into specific radial ranges marked by the three anchor points 
for each galaxjQ. The resulting color profile slope represents 
the change in color (A mag) per dex in radius, where a positive 
slope indicates that the galaxy is getting redder with increasing 
radius from its center 

We implement a Monte Carlo method to compute the color 
gradients and their uncertainties. The essence of our ap- 
proach is to iteratively sample the two observables, the radius 
and the color at that radius, while incorporating uncertainties 
from Poisson noise, sky subtraction, radius measurement, and 
stochastic fluctuations in the color profile (due to, for instance, 
substructure from dust lanes, star clusters, or spiral arms). Ex- 
cept for Poisson noise, which is random, the other uncertainties 
generally contribute to errors in a systematic way. Our Monte 
Carlo approach makes it possible to derive an effective random 
uncertainty from a series of systematic errors. We start with the 
premise that the effective radius is uncertain by a Gaussian er- 
ror with (J w QAR^o, which is frequently true when the data are 
not in the noise-dominated regime. The systematic uncertainty 
on /?5o arises from the fact that there is a range of plausible, 
acceptable models that scatter around the best-fitting solution. 
We draw a radius from that distribution, and the color at that 
corresponding radius. The color value is also sampled from a 
Gaussian distribution given by the Poisson noise, centered on 
the original color at the sampling radius. The net random un- 
certainty of the color includes uncertainty in the sky value. The 
same sampling process is applied at l.SRso, now with an ef- 
fective scatter of l.SaiR^o), and in the central region with a 
scatter of the seeing FWHM. Having sampled the radii and col- 
ors around Rq, R^o, and 2.5/?5o, we calculate the color gradient 
following 



V(Color) 



in/out ■ 



(Color),, -(Color),, 



(7) 



logri-logro 

This procedure is repeated lO"* times to generate a distribution 
of color, the median and width of which are the gradient and 
the uncertainty in the gradient, respectively. 

Figure |6] shows the normalized histograms of the inner and 
outer gradients of the B-I color, sorted into bins of differ- 
ent morphological types; the corresponding gradients for the 
V-I and R-I colors are shown in Figures |7]and|6] respectively. 
We define positive gradients as those larger than 0. 1, flat gradi- 
ents as those between -0.1 and 0.1, and negative ones as those 
smaller than -0. 1 . The vertical dashed lines in each panel mark 
the boundaries between these three categories. 

Table 3 summarizes the statistics for CGSFl The color gradi- 
ents depend strongly on galaxy morphology and color. Ellipti- 
cal galaxies generally have very little, if any, measurable color 
gradients. Interestingly, their central regions show a slight ten- 



dency to exhibit positive gradients in all three colors, whereas 
beyond their effective radii the trend reverses and there is a mild 
preference for negative gradients. In either case, the distribution 
of gradients is narrowly peaked, with a dispersion of ^ 0. 1 1 . SO 
and so/a galaxies largely follow the same pattern as the ellipti- 
cals. By contrast, spirals of types Sa through Sd behave quite 
differently. The inner regions of these galaxies show a wide 
dispersion in gradients (^ 0.23), and they are predominantly 
negative: the colors get redder toward the center. The gradients 
in the outer regions, on the other hand, are predominantly flat 
(peak near 0), and there is roughly an equal number of positive 
and negative values, although the scatter is large. Galaxies be- 
longing to the latest types (Sdm and Sm) display no preference 
for gradients of either sign, neither in their interior nor in their 
exterior regions. The above trends stand out most clearly in 
B-I, the color combination with the greatest wavelength sep- 
aration, and they become less pronounced in V-I, and even 
more so in R-I, although they are still noticeable. 

7. BARS 

Two fundamental quantities that characterize a bar are its 
length and strength. There are many ways to estimate the 
characteri stic size of a bar . Apart from simple visual in- 
spection (iKormendvl Il979l) . the most commonly used ap- 
proaches involv e measurement of the maximum value of the bar 
ellipti city (e.g., iLaine et akl 120021: iMenendez-Delmestre et al 



2005 



2007 ), the radial variation of the positi on angle (e.g.. lErwin 



iMenendez-Deknestre et al. I l2007h . the radial variation 

of the phase angle of the second Fourier mode (Agu erri et al.l 



[2003), and the bar-interbar contrast (Aguerri et ak 12003^ 
as well as direct decomposi t ion o f the image into differ- 
ent components (iPrieto et al.l Il997l) . The strength of the 
bar can be ascertaine d by quantifying the maximum elliptic- 
ity in the bar r egion (Martin 'l995t [Martinet & Friedlil [T997t 
IMenendez-Delmestre et al. 2007 ). the torques generated by the 
bar ( Buta & Blockl 200 ll) . the amplitude of the even Fourier 
modes of the isophotal intensity distribution ( lOhta et al.ll 19901: 
lAthanassoula & Misiriotisll2002l) . and direct decomposition of 
the galaxy into its co nstituent light fra ctions (iLaurikainen et al.l 
120051: lGadottill2008l: iPeng et al.ll2010h . Here we present a pre- 
liminary appraisal of the bar properties of the CGS sample 
based on information that can be readily extracted from our 1- 
D isophotal data. We describe analyses based on the geometric 
parameters and Fourier components. 

7.1. Geometric Analysis 

In the absence of confusion from dust, star-forming regions, 
and projection effects, bars usually leave a distinctive imprint 
on the e and PA profile of a galaxy. The bar is marked by a re- 
gion wherein the ellipticity rises steadily until it reaches a peak 
and drops, and, unless the bar semi-major axis is closely aligned 
with the major axis of the disk, the constant position angle in 
the bar regi on abruptly changes value as it transitions into the 
disk region (iGadotti et al.ll2007l) . 

We begin with the e and PA profiles of the /-band image, 
as extracted from the second step of running the task ellipse 
(Section 4.1), during which only the galaxy center was held 
fixed. The / band is preferred over the other bands at shorter 
wavelengths because it mitigates contamination by dust and 
young stars. We consider the profiles from an inner radius 



'in the future, we will also derive the color gradient in physically interesting regions, such as spiral arms, bars, or the break points of the surface brightness profiles. 
'"We excluded the galaxies flagged in Table 2 as adversely affected by field stars. 
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TABLE 3: Frequency of Color Gradients 



Hubble Type 


T 




Color 




Inner Region 






Outer Region 














Positive 


Flat 


Negative 


Positive 


Flat 


Negative 


(1) 


(2) 




(3) 




(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


E 


T < 


-3 


B - 


I 


12.0% 


85.5% 


2.5% 


4.8% 


36.1% 


59.1% 








V - 


I 


58.0% 


42.0% 


0.0% 


13.6% 


50.6% 


35.8% 








R- 


I 


65.9% 


34.1% 


0.0% 


12.2% 


65.9% 


21.9% 


SO-SO/a 


— 3 < T 


< 0.5 


B - 


I 


16.2% 


63.5% 


20.3% 


27.0% 


24.3% 


48.7% 








V - 


I 


35.1% 


57.1% 


7.8% 


31.2% 


33.8% 


35.0% 








R- 


I 


49.3% 


50.7% 


0.0% 


32.0% 


49.3% 


18.7% 


Sa-Sab 


0.5 < T 


< 2.5 


JD — 


T 
1 


7.8% 


27.5% 


64.7% 


15.7% 


23.5% 


60.8% 








V - 


I 


9.8% 


58.8% 


31.4% 


23.5% 


25.5% 


51.0% 








R- 


I 


17.6% 


78.4% 


4.0% 


23.5% 


52.9% 


23.6% 


Sb-Sbc 


2.5 < T 


< 4.5 


B - 


I 


1.5% 


10.7% 


87.8% 


35.9% 


10.7% 


53.4% 








V - 


I 


2.3% 


39.1% 


58.6% 


27.8% 


24.1% 


48.1% 








R- 


I 


3.0% 


85.8% 


11.2% 


40.3% 


33.6% 


26.1% 


Sc-Sd 


4.5 < T 


< 7.5 


B - 


I 


6.0% 


8.2% 


85.8% 


47.0% 


10.4% 


42.6% 








V - 


I 


4.4% 


37.0% 


58.6% 


43.0% 


20.0% 


37.0% 








R- 


I 


3.7% 


80.6% 


15.7% 


50.0% 


26.1% 


23.9% 


Sdm-Sm 


7.5 < T 


< 9.5 


B - 


I 


23.5% 


41.2% 


35.3% 


58.8% 


17.6% 


23.6% 








V - 


I 


17.6% 


70.6% 


11.8% 


41.2% 


29.4% 


29.4% 








R- 


I 


12.5% 


81.3% 


6.2% 


50.0% 


31.3% 


18.7% 



NOTE. — Column 1: Hubble type. Column 2: morphological type index. Column 3: color. Columns 4-9: fraction of the 
positive, flat, and negative color gradients in the inner and outer regions. 



corresponding to 3 times the seeing disk to the radius where 
the isophotal intensity is 1 cr above the sky background, be- 
y ond which we truncate the surface brig htness profile. Simila r 
to lMenendez-Delmestre et"ar ] ( l2007h and lAguerri etal] (l2009l) . 
bars are required to have a maximum projected ellipticity (cmax) 
greater than 0.2, and within the bar region the position angle 
should be constant to within APA < 20°. If none of the data 
points in the e profile exceeds 0.2, or if Ae < 0.1 through- 
out the entire e profile, we classify the galaxy as unbarred. 
If Ae > 0.1 somewhere along the e profile but the associated 
APA < 10°, it is possible that a bar exists but happens to 
align fortuitously with the major axis of the outer disk. We 
flag these cases as "possibly" barred and carefully inspect the 
galaxy image visually to see if we can confirm their reality. If 
the galaxy is barred, we set the inner boundary of the bar re- 
gion to be the first data point in t he e profile that exceeds 0.2. 
iMenendez-Delmestre et al.l (l2007h find that near the end of the 
bar the ellipticity and position angle usually begin to show large 
deviations, typically at the level of Ae > 0.1 and APA > 10°. 
We adopt these criteria to define the radius of the outer bound- 
ary of the bar. 

The projected bar size is set to be the semi-major axis o f 
the isophote where e peaks ( IMenendez-Delmestre et aL I I20071) . 
with the associated error as the semi-major axis range contain- 
ing the tip of the e profile (i.e. e > emax-O.Ol) in the bar re- 
gion. The reason we do not simply equate the bar size with the 
outer boundary of the bar is because the change in e and PA 
in the transition zone between the bar and spiral arms can be 
influenced by the latter The outer boundary of the bar can be 
overestimated if the bar is aligned with the spiral arms. This 
effect can be mitigated by using the position where e peaks as 
the bar size, since spkal features cannot produce ellipticity val- 



ues as high as those of a bar. Assuming that the intrinsic shape 
of the galaxy disk is purely circular, we correct the observed, 
projected bar length (R^,^^) to its intrinsic value (R\,^^) following 

„ I ' TTsinAPAV 
^bar = ^Ly(cosAPA)2+(^^— — j , (8) 

where APA = PAg^i - PAbai , and PAgai and egai are the position 
angle and ellipticity of the outer disk of the galaxy, respec- 
tively. Although a full 2-D analytical deprojection of the bar 
is more accurate (see Appendix A in iGadotti et al.ll2007h . in 
practice the size estimates from the two methods agree very 
well. The typical difference in bar radii measured by 1-D fit- 
ting compared to 2-D deprojection is about 0"3. Since the 
difference is very small, we use the 1-D method for simplic- 
ity. The projected emax then repre sents the strength of the bar 
( IMenendez-Delmestre et alj|2007h . with the fitted error of e^ax 
as its uncertainty. The position angle of the bar, PAbai , is given 
by the average PA over the bar region, with the rms as its error. 

Figure|9|illustrates how we identify and measure the bar size 
and strength, as applied to the star-cleaned /-band image of the 
SBb galaxy NGC 7513. The radial profiles of e and PA clearly 
show the hallmark features of a bar: a distinctly broad peak in e 
above our minimum threshold of 0.2, reaching e^ax = 0.68, and 
an extended plateau of near-constant PA « 75° (APA < 20°). 
The outer disk of the galaxy has a clearly different e (0.33) 
and PA (105°). The two vertical solid lines mark the inner 
and outer boundaries of the bar-dominated region; the verti- 
cal dotted line gives the projected bar radius. Although widely 
used in the literature, the measured e^ax of the bar is actually 
~ 20% lower than that derived from 2-D image decomposition 





Fig. 9. — Illustration of how we determine the size and strength of a bar. Left: star-cleaned /-band image of NGC 7513; the size of the image is ~3' x3'. Right: 
radial profiles of PA, e, h/h, and 4'2- The horizontal dashed lines in the PA and e panels denote the characteristic values of the galaxy. The solid vertical lines 
mark the inner and outer boundaries of the bar-dominated region, and the dotted vertical line represents the projected size of the bar determined using the geometric 
method, where the bar is required to have emax > 0.2 and APA < 20°. In the /2//0 and (/>2 panels, the soHd vertical line marks the inner boundary of the bar. 
The vertical dashed line represents both the outer boundary and size of the bar, determined using the Fourier method, where the bar criteria are h/Io > 0.2 and 
A02 < 20° . The corresponding isophotal ellipses for the two methods are overplotted on the left-hand image with the same type of lines as in the right-hand panel. 
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Fig. 10. — Left: correlation between R^.J'^'"'', the deprojected bar size measured from the geometric, ellipse method and /?Jj™"", that measured from the Fourier 
decomposition. The dashed line represents y = x. Right: correlation between </>2 and PA of the bar. The dashed line represents y = 180° - .v. 



when the bulge and disk components are also included jGadottil 
l2008h . Indeed, reliable bar parameters can only be determined 
by 2-D decomposition of the images with all of the other com- 
ponents included. Such analysis is outside the scope of this 
work, but future papers will present results from 2-D decompo- 
sitions of the CGS images. 

7.2. Fourier Analysis 

In addition to the geometric method described above, we 
also derive bar properties using the radial profiles of the rela- 
tive amplitude of the m = 2 Fourier mode (h/h) and its asso- 
ciated phase angle (</)2). As before, we work with the /-band 



images. Bars are usually associated with the first local maxi- 
mum in the h/h profile, where the bar/interbar contrast is the 
largest, over an extended region where 02 keeps approximately 
constant. Subsequent maxima in the h/h profile, if present, 
trace spiral arms or ring structures, but in these instances 4>2 
varies with radius. Spiral arms always produce varying phase 
angles, and thus the region where they dominate can be easily 

excluded from the bar size measurement. 

Ado pting a procedure similar to that used by lAguerri et al.l 
(1200^, we define the bar to be the region wherein the maxi- 
mum relative m = 2 Fourier amplitude (h/Iojm-dx > 0.2 and the 
phase angle remains constant to A(l>2 < 20°. We set the inner 
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TABLE 4: Properties of Bars 



Name 


Flag 




Cbar 


PAbar 


(-^2/Jo)bar 


(02) bar 


Notes 






(") 




n 




(°) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


ESO 009-GOlO 


N 














ESQ 027-GOOl 


B 


23.74±4.45 


0.50±0.02 


64.25±3.77 


0.38±0.04 


126.77±6.04 




ESO 027-G008 


N 












Extra 


ESO 056-G115 


N 












Bvnr 


ESO 060-G019^ 


B 


20.43±3.00 


0.38±0.02 


28.58±5.00 


0.76±0.03 


145.03±5.56 




ESO 091-G003 


N 












Extra 


ESO 097-G013 


N 












BVRr 


ESO 121-G006 


N 














ESO 121-G026 


B 


19.21±2.26 


0.57±0.01 


166.75±11.08 


0.55±0.03 


176.96±1.45 




ESO 136-G012 


N 












BVRr 



NOTE. — Column 1: galaxy name. Column 2: flag indicating whether the galaxy is barred: "B" = barred; "B:" = possibly 
barred; "N" = not barred. Column 3: bar radius (intrinsic, deprojected according to Equation (7)). Column 4: bar ellipticity. 
Column 5: position angle of the bar, east of north. Column 6: relative amplitude of the m = 2 mode of the bar. Column 7: phase 
angle of the m = 2 mode of the bar. Column 8: notes: "extra" = object not formally part of the main sample; filters = image of 
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Fig. 11. — Illustration of how we measure the strength of the spiral arms using the Fourier method. Left: star-cleaned /-band image of NGC 5247; the size of the 
image is ~7f3 x7f3. Right: radial profiles of PA, e, /2//0, and ipi- The solid vertical lines mark the inner and outer boundaries of the disk-dominated region, which 
is defined to be outside the central bulge or bar component, but inside the radius where the isophotal intensity is 3 cr above the sky background. The corresponding 
isophotal ellipses are overplotted on the left-hand image. The horizontal dashed lines in the PA and e panels denote the characteristic values of the galaxy. 



boundary of the bar to be the first data point outside of 3 times 
the seeing radius in which h/h) > 0.2. Past the peak, we 

designate the radius where h/Io = (/2/A))max-0.1 as the outer 
boundary of the bar region. In the event that there is a secondary 
maximum and the local minimum between the two peaks ex- 
ceeds (/2//())max-0.1, we sct the position of the local minimum 
to be the outer boundary of the bar region. 

We assign (/2//o)max to be the bar strength, with the uncer- 
tainty set by the statistical error derived from the Fourier de- 



composition process. As the Fourier method is minimally af- 
fected by spiral arms, we simply set the bar size to be equal to 
the radius of its outer boundary; its associated uncertainty is the 
semi-major axis range between the outer boundary and the ra- 
dius where h/Io = (h/Io) max -0.05. The Fourier analysis is per- 
formed on the isophotes extracted in the third step of the ellipse 
fitting, where the geometric parameters were all fixed to those 
of the outermost isophote (see Section 4.1). For simplicity, we 
just assume that the disk is purely circular in its face-on view. 



Fig. 12. — Illustration of how we measure the strength of the spiral arms using the structure map. Left: star-cleaned composite color image of NGC 5247; the 
size of the image is ~7'.3x7'.3. Right: structure map of the star-cleaned B-band image. The two overplotted ellipses mark the inner and outer boundaries where the 
spiral arms populate. After rejecting the masked objects, we use the rms of the pixels within these two ellipses to estimate the strength of the spiral arms. (A color 
version of this figure is available in the online journal.) 



SO the semi-major axes of the isophotal ellipses are the radii of 
the circles. The size of the bar is denoted by the semi-major 

axis of the particular isophote that encloses the bar region, 
which is actually the radius of the circle that passes right 
through the end of the bar in a face-on view. Under this sim- 
plified assumption, the bar size is already deprojected. The 
characteristic phase angle of the bar is the average 02 within 
the bar region, with A(f>2 as the uncertainty. The two bottom- 
right panels of Figure |9] show our Fourier technique applied to 
NGC 7513. 

Fourier analysis offers a very useful approach to studying 
bars. It not only provides another independent method to iden- 
tify and quantify bars, but also recovers bars missed by the 
geometric method in galaxies with internal structure too com- 
plex to yield an unambiguous bar signature in their e and PA 
profiles. In fact, the bar parameters derived from these two 
methods usually agree quite well. Figure [TO] (left panel) com- 
pares the bar sizes measured from the ellipse method (TJ^^J.'''*'''; 
deprojected) versus those measured from the Fourier method 
(^bar""''')- Overall, there is a good correlation, but some glar- 
ing outliers stand out. Although our analysis is done in the / 
band, dust extinction can still be significant in some galaxies. 
The dust extinction features near the central region of a galaxy 
tends to trick the Fourier method by producing a false peak in 
h/Io with roughly constant this yields an unusually com- 
pact, incorrect bar size. The Fourier method can also be unreli- 
able for weak bars embedded in disks with strong spiral arms. 
Under these circumstances, the relatively weak local peak of 
the bar in the I2 / lo profile can be overshadowed by a stronger 
peak generated by the spiral arms, leading to an overestimate 
of the bar size. The position angle of the bar correlates strongly 
with 02 (Figure[TOl right panel). Ideally, 02 = 180° -PA. How- 
ever, when the PA of the bar approaches 0° or 180°, the cor- 
responding phase angle can have values similar to the PA; this 



is responsible for the points lying on the lower-left and upper- 
right portions of the figure. Dust in the central regions of the 
galaxy further contributes to the scatter 

7.3. Final Bar Classification 

Since neither of the methods discussed above is foolproof, 
we use both to assign the final bar classification to the galaxies 
in CGS. If a galaxy is not classified as barred in both the geo- 
metric and Fourier analysis, it is labeled unbarred. If the galaxy 
is classified as barred by only one method, we call it "possibly" 
barred. If both methods consider a galaxy barred, we check 
whether the derived bar sizes are consistent between them. We 
only classify it as definitely barred if the bar sizes from the two 
methods differ by less than 10"; if they disagree by more than 
10", we call it "possibly" barred. 

For the "possibly" barred galaxies, we visually examine 
their /-band images and further scrutinize their geometric and 
Fourier profiles to check whether we are being misled by inter- 
nal structural complexities such as spiral arms or dust features. 
If any of the measurements from either of the two methods is 
suspect, we manually set the inner and outer boundaries of the 
bar region and redo the measurements. Not all ambiguous cases 
can be resolved, and in our final classification we continue to 
flag their bar status as uncertain. For our final measurement of 
the bar size, we usually give higher priority to the results de- 
rived from the geometric method, which generally gives more 
accurate sizes than those derived from the Fourier method. The 
definition of the bar size in the Fourier method is of limited 
utility because it is largely arbitrary and more prone to being 
influenced by the presence of spiral arms. Table 4 summarizes 
the final bar classification and some basic bar parameters de- 
rived from our analysis. Among the 501 disk galaxies {T > -3) 
in the final catalog, 44 (9%) are deemed definitely barred, 136 
as possibly barred (27%), and 321 (64%) as unbarred. As bar 
identification is uncertain in highly inclined galaxies, we reex- 
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amine the statistics in a subsample restricted to have ellipticities 
smaller than ggai = 0.6. As expected, the bar fraction increases. 
Out of 387 disk galaxies, 173 (45%) are barred and 214 (55%) 
are unbarred. A more detailed comparison between our results 
and those in the literature will be deferred to a future paper. 



a large I\/k), but 0i increases monotonically as a function of ra- 
dius. 

Our method to measure lopsidedness differs somewhat from 
that used in previous works. Rix & Zaritsky ( 1995j) perform a 
bulge-to-disk decomposition of the surface brightness profile to 
determine the scale length of the disk, and then compute the 
average relative Fourier amplitude between 1.5 and 2.5 scale 
lengths. As we do not yet have robust structural decomposi- 
tions for the entire CGS sample (this work is in progress), we 
resort to a simpler strategy, one based on the expectation that 
the lopsided portion of the disk should be characterized by a 
roughly flat 01 radial profile. Through careful experimentation 
with a number of galaxies with prominent lopsided disks, we 
find that the lopsided region can be effectively isolated by re- 
quiring that the phase angle be constant to within A^i < 70°. 
We apply this criterion to the phase angle profile of the /-band 
image to define the inner and outer radii of the lopsided region, 
and then adopt these values for the other filters. For each filter, 
the lopsidedness is the average value of I\ //q within that region, 
with the standard deviation as its associated uncertainty. The 
characteristic value of 0i and its error are calculated similarly. 
Figure[T4]illustrates our method on the Scd galaxy NGC 7070. 

Lopsidedness measurements for the entire sample are given 
in Table 43- Table 7 summarizes the frequency of galaxies 
with significant lopsi dedness, defined as (/i//o), > 0.2 (e.g., 
IZaritskv & Rixl [19971) . for the subsample of 350 disk galaxies 
with ggai < 0.6 for which robust measurements could be made. 
Consistent with previous studies, the fraction of galaxies with 
significant lopsidedness is high, and the frequency increases 
substantially in late-type galaxies. 



8. SPIRAL ARMS 

We provide two quantitative measurements (Table 5) that can 
be used to access the presence and strength of spiral arms in 
galaxies. The analysis is applied uniformly to all non-elliptical 
galaxies in the sample, including disk galaxies traditionally 
deemed to lack spiral arms, such as SOs. 

We perform a simple measurement of the average strength of 
h/k) in the disk region outside the central bulge and the bar. 
Spiral arms are the main contributor to any significant m = 2 
mode in this regiorQ If neither a featureless, classical bulge 
nor a bar is present, the minimum inner boundary for our cal- 
culation is set to 3 times the seeing radius. For barred galaxies, 
the inner boundary is naturally set to the bar radius, which al- 
most always lies exterior to the bulge. For unbarred galaxies 
with classical bulges, we define the inner boundary to be the 
radius where e > 0.2, beyond which the disk usually dominates 
over the bulge. This criterion fails for face-on galaxies with 
very weak spiral arms and classical bulges, because the disk 
becomes indistinguishable from the bulge on the basis of its el- 
lipticity alone. Fortunately, under these circumstances both the 
bulge and the disk contribute little to /2//0 anyway, and it makes 
little difference whether the bulge is excluded or not. The outer 
boundary is the radius where the isophotal intensity reaches 3 cr 
above the background in the /-band image; we apply the same 
boundary for the other filters. We then calculate a characteris- 
tic value of I2 //() by averaging its profile between the inner and 
outer boundaries. We illustrate our methodology in Figure [TT] 
applied to the Sbc galaxy NGC 5247. 

Our second method makes use of the structure maps (Paper I) 
to estimate the strength of the spiral features (Figure[T2li. After 
masking out the field stars and background galaxies, we com- 
pute the standard deviation of all the remaining pixels within 
the inner and outer boundaries of the disk-dominated region, 
as determined above. The agreement between the two different 
measurements is not good, as can be seen in Figure [T3] where 
we plot (/2//0) against as for all the filters. The two parame- 
ters trace structures on different scales. The structure map op- 
timally filters spatial features on the scale of the point-spread 
function. It effectively highlights features such as dust lanes 
and thin arms, but it is not very sensitive to smooth and wide 
spiral arms, which can be better probed via {Ii/h). 

9. LOPSIDEDNESS 

The relative amplitude of the m = 1 Fourier mode is 
widely used to study the lopsidedness of galactic stellar 
jRix & Z aritsky '1995; 'Zaritskv & Rix' '1997'; 'Bo urnaud et all 
I2OO5I; iJog & Combes 200 8; Re ichard et a l. 2008) and gaseous 
(e.g., Ivan Evmeren et al1l201 ll) disks. This approach is well- 
defined, quantitative, and relatively straightforward to imple- 
ment, and hence can be applied to study large samples of galax- 
ies. A lopsided disk stands out as a region of enhanced /i //q and 
roughly constant phase angle <f>i . A one-arm spiral also exhibits 

"The m = 2 mode is most sensitive to systems with grand design, two-arm spirals. However, flocculent or multiple-arm spirals still exhibit significant amplitude in 
the m = 2 mode. We defer a full treatment of spiral arms, including exploration of high-order modes, to a separate paper 

'^We have flagged the galaxies whose lopsidedness measurements may be unreliable because of contamination by bright stars or excessive crowding by field stars. 



10. DATA VERIFICATION 
10.1. Internal Comparison 

Several galaxies were observed more than once during differ- 
ent nights throughout the survey. Although only the best images 
are included in the final CGS catalog, the repeat observations, 
which were reduced and analyzed in the same manner as the 
rest, afford an opportunity to access the accuracy of our calibra- 
tion methods and the reliability of the parameter measurements. 

Table 8 lists the galaxies that have pairs of repeat observa- 
tions useful for internal comparison. Note that for this exercise 
we only select objects that have reasonably good data. Many 
of the duplicate observations were taken precisely because the 
original observation was deemed to be of exceptionally low 
quality, either because of weather conditions (bad seeing, ex- 
cessive cloud cover) or technical problems (poor telescope fo- 
cus, tracking errors). Because one of the observations in the 
comparison pair is — by definition — suboptimal, the following 
assessment, in some sense, gives an overly conservative esti- 
mate of the magnitude of internal errors. 

Figure[T5]compares the surface brightness profiles for 
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TABLE 5: Properties of Spiral Features 



Name 


{I2/I»)b 














4 


Notes 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


ESO 009-GOlO 


0.11±0.08 


0.09±0.07 


0.08±0.06 


0.07±0.05 


0.006 


0.004 


0.003 


0.002 




ESQ 027-GOOl 


0.31±0.21 


0.29±0.20 


0.28±0.21 


0.27±0.23 


0.008 


0.008 


0.005 


0.002 




ESO 027-G008 


0.23±0.19 


0.22±0.18 


0.23±0.18 


0.23±0.20 


0.010 


0.008 


0.007 


0.003 


Extra 


ESO 056-G115 


















BVRI^ 


ESO 060-G019 


0.55±0.24 


0.50±0.23 


0.48±0.21 


0.48±0.16 


0.044 


0.038 


0.048 


0.037 




ESO 091-G003 


0.10±0.06 


0.11±0.06 


0.11±0.07 


0.12±0.08 


0.008 


0.008 


0.005 


0.003 


Extra 


ESO 097-G013 


















BVRI^ 


ESO 121-G006 


0.25±0.22 


0.24±0.20 


0.22±0.17 


0.18±0.10 


0.009 


0.020 


0.027 


0.025 




ESO 121-G026 


0.20±0.09 


0.19±0.09 


0.19±0.08 


0.17±0.07 


0.010 


0.008 


0.007 


0.003 




ESO 136-G012 


0.22±0.07 


0.33±0.58 


0.23±0.28 


0.19±0.15 


0.013 


0.017 


0.016 


0.008 


BVRI'° 



NOTE. — Column 1: galaxy name. Columns 2-5: average value of /2//0 within the disk-dominated region of the B— ,V — , R— , 
and 7-band images, respectively. Columns 6-9: standard deviation of the pixels within the disk-dominated region of the structure 
maps of the B—,V—,R—, and 7-band images, respectively. Column 10: notes: "extra" = object not formally part of the main 
sample; filters = image of that filter has problem, which is indicated by the table mark. (This table is available in its entirety in a 
machine-readable form in the online journal. A portion is shown here for guidance regarding its form and content.) 

"Galaxy adversely affected by a very bright star. 

''Galaxy adversely affected by very crowded field stars. 

■^Galaxy distorted by an interacting neighbor. 



TABLE 7: Frequency of Significant Lopsidedness 



Hubble Type 


T 


Number 


Fraetiuu 


(1) 


(2) 


(3) 


(4) 


so-so/a 


-3<T<0.5 


77 


15.6% 


Sa-Sab 


0.5<T<2.5 


48 


25.0% 


Sb-Sbc 


2.5 < T < 4.5 


110 


45.5% 


Sc-Sd 


4.5 < T < 7.5 


103 


64.1% 


Sdm-Sm 


7.5 < T < 8.5 


12 


83.3% 



Note. — Cohimu 1: Hubble type. Column 2: morphological type index. Colrmur 
3; number of objects. Colurrm 4; fraction of galaxies with significant lopsidedness, 
defined to be those with {h/lo) i > 0.2 in the 1 band. 



galaxies observed on different nights. In most cases, they 
agree quite well. The weighted average of the profile dif- 
ferences (dashed line, calculated by averaging the profile dif- 
ferences weighted by the corresponding error at each data 
point) resides well inside the formal 1 a photometric uncer- 
tainty (iNoordermeer & van der Hulsjl2007h . This suggests that 
our photometric errors are robust, both for the photometric 
and non-photometric observations. A few objects (NGC 1374, 
2196, 6810, 6861, 7590) show slightly larger, but by no means 
alarming, discrepancies. Two types of differences in profile 
shape can be seen. The innermost portions of the profiles of- 
ten show systematic deviations, sometimes as large as ^ 0.5 
mag arcsec"^. This effect can be entirely attributed to mis- 
matches in seeing, but is well confined within ~3 times the 
radius of the seeing disk. This is the reason we restrict all of 
our scientific analysis to radii beyond this. Additionally, many 
of the profiles show some level of systematic deviation at large 
radii. This most likely arises from errors in sky subtraction. 
For most objects, the deviations occur at the level of ~ 0.2 
mag arcsec"^, but they lie well within the error bars of the indi- 
vidual isophotal intensities, which again indicates that our error 
budget is realistic. The most extreme deviations occur in galax- 
ies that are too extended 



TABLE 8: Internal Comparison Sample 



Name 


Filter 


Date 1 


Conditions 1 


Date 2 


Conditions 2 






(yyyy-mm-dd) 




(yyvv-mm-dd) 




(1) 


(2) 


(3) 


(4) 


■ ■ ■ ' (5) 


(6) 


ESQ 311-G012 


BVRI 


2006-01-28 


N 


2004-01-18 


N 


IC 2597 


BVRI 


2006-01-27 


N 


2004-04-16 


N 


NOO 1374 


VI 


2003-11-20 


N 


2003-11-22 


P 


NOO 1399 


BVRI 


2006-01-27 


N 


2003-11-18 


N 


NGC 1688 


BVRI 


2006-01-28 


N 


2004-01-17 


N 


NGC 1700 


BVRI 


2006-01-27 


N 


2003-11-20 


N 


NGC 2196 


BVRI 


2006-01-29 


N 


2003-02-09 


N 


NGC 2310 


BVRI 


2006-01-29 


N 


2004-01-18 


N 


NGC 2434 


BVRI 


2006-01-26 


N 


2003-11-20 


N 


NGC 2698 


BVRI 


2006-01-29 


N 


2004-04-16 


N 


NGC 3672 


HI 


2004-01-20 


P 


2004-01-17 


N 


NGC 6810 


V 


2003-05-28 


P 


2005-09-04 


N 


NGC 6861 


BVRJ 


2005-09-08 


N 


2003-05-25 


N 


NGC 6887 


V 


2005-05-07 


P 


2005-09-06 


N 


NGC 6943 


BVRI 


2005-09-08 


N 


2005-05-10 


N 


NGC 7083 


BVRI 


2005-09-08 


N 


2004-10-14 


N 


NGC 7329 


BVRI 


2005-09-05 


N 


2004-10-14 


N 


NGC 7410 


V 


2005-09-06 


N 


2005-09-09 


N 


NGC 7590 


BVRI 


2005-09-08 


N 


2004-10-10 


N 



NOTE. — Column 1: galaxy name. Column 2: filter. Columns 3 and 4: observation date and 
photometric condition of tiie images used in the primary sample. Columns 5 and 6: observation 
date and photometric condition of the images used in the comparison sample. In Cols. (4) and 
(6), "P" means photometric, and "N" means non-photometric. 



for standard sky subtraction, for which we had to resort to an 
indirect estimate based on profile fitting (Section 3). In these 
cases, the deviations in the outer profiles may be as large as 
^ 0.5 mag arcsec"^. NGC 2434 is such an example. However, 
our formal error bars appear to be realistic even in these extreme 
situations. 

Figure [16] compares nine measured parameters derived from 
the set of repeat observations. Observation 1 denotes the mea- 
surement with better quality that has been adopted in the final 
database of the survey, and Observation 2 gives the compar- 
ison measurement. We can see that overall the agreement is 
quite good. Notable exceptions can be identified with galaxies 
that have especially unreliable sky values, such as NGC 2434, 
which is the most deviant outlier in the R^q plot (this quantity is 
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0.01 



0.02 



0.03 







0.01 0.02 0.03 0.04 



Fig . 13. — Comparison between the two measures of the spiral arms ( (h / If, ) and CTs ) for BVRI filters. The correlation is very weak or absent, because the two 
measures essentially probe structures on quite different scales. CTs is mainly determined by structures with typical scales of the point-spread function, while {h/kt) 
is more sensitive to more extended, large-scale features coherent over significant portions of the spiral aims. 



sensitive to data at large radii). The average differences and 
standard deviations of the parameters plotted in Figure [16] are 
listed in Table 9. The average differences are quite close to 0, 
and the scatter is small. 

10.2. External Comparison with HyperLeda 

Several of the global parameters measured i n our sample 
have independent data listed in HyperLed^H dPaturel et al.l 
I2003h . which we can use to perform an external comparison 
of our errors. Figure[T7]compares the following six parameters 
between CGS and HyperLeda: total /-band magnitude (/tot), to- 

^ ^http : / / leda . univ-lyonl . f r| 



tal B-band magnitude (/?tot), integrated B-V color, isophotal 
diameter at 25 B mag arcsec"^ {D25), e, and PA. The overall 
agreement is quite good for the integrated magnitudes, B-V 
color, and D25 diameters, especially for the red open points, 
which represent galaxies that were observed under photometric 
conditions, that have no contamination from the nearby bright 
stars, and that have reliable sky values. The PA comparison 
improves dramatically after isolating the subset with e > 0.3. 
When the galaxy is round, the PA is hard to determine because 
the semi-major axis for any given isophote is ill-defined. In ad- 
dition, extreme outliers lying on the lower-right and upper-left 



Fig. 14. — Illustration of how we measure lopsidedness. Left: star-cleaned /-band image of NGC 7070; the size of the image is ~3f6x3f6. Right: radial 
profiles of PA, e, I\/Iq, and 0i . The horizontal dashed lines in the PA and e panels denote the characteristic values of the galaxy. The solid vertical lines, and the 
con'esponding isophotal ellipses oveiplotted in the left-hand image, mark the inner and outer boundaries of the region used to compute the lopsidedness, which we 
define to be that in which the radial vaiiation of 0i is smaller than 70°. The lopsidedness is measured simply by averaging /[ //q within this region. 



TABLE 6: Lopsidedness Measurements 



Name 


{4'i)i 




{h/h)y 






Notes 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


ESO 009-GOlO 


338.5±9.9 


0.20±0.10 


0.16±0.09 


0.16±0.09 


0.15±0.10 




ESQ 027-GOOl 


70.4±4.0 


0.53±0.14 


0.50±0.14 


0.48±0.12 


0.45±0.11 




ESO 027-G008 


168.9±2.3 


0.13±0.03 


0.10±0.03 


0.14±0.05 


0.21±0.13 


Extra 


ESO 056-G115 














ESO 060-G019 


240.0±11.6 


0.12±0.10 


0.18±0.08 


0.16±0.05 


0.19±0.05 




ESO 091-G003 


11.0±17.4 


0.10±0.04 


0.09±0.02 


0.10±0.04 


0.17±0.09 


Extra 


ESO 097-G013 












BVRf' 


ESO 121-G006 


226.3±18.4 


0.25±0.09 


0.29±0.12 


0.28±0.07 


0.17±0.06 




ESO 121-G026 


161.0±19.6 


0.30±0.13 


0.35±0.15 


0.34±0.15 


0.29±0.14 




ESO 136-G012 


155.3±17.0 


0.23±0.19 


1.01±1.62 


0.42±0.58 


0.28±0.25 


BVRI'° 



NOTE. — Column 1: galaxy name. Column 2: phase angle (in units of degrees) of the lopsided region in the /-band image. 
Columns 3-6: strength of the lopsidedness parameter for the B—,V—,R—, and /-band images, respectively. Column 7: notes: 
"extra" = object not formally part of the main sample; filters = image of that filter has problem, which is indicated by the table 
mark. (This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance 
regarding its form and content.) 

"Galaxy adversely afl'ected by a very bright star. 

''Galaxy adversely affected by very crowded field stars. 
Galaxy distorted by aji interacting neighbor. 



corners of the distribution can be attributed to the 180° ambi- 
guity for PAs close to 0° or 180°. The ellipticities show by far 
the worst agreement. The large scatter may be due in part to the 
fact that the HyperLeda values pertain to measurements made at 
Hb = 25 mag arcsec"^, whereas ours are made at a significantly 
lower surface brightness threshold, at /ig « 27 mag arcsec"^. 
A more serious problem may be related to inherent biases and 
errors that are known to plague the axial ratios (ellipticities) 
contained in the HyperLeda database (see the Appendix of Ho 



2007). 

10.3. External Comparison with SDSS 

Roughly 9% of the CGS galaxies overlap with the Sloan 
Digital Sky Survey (SDSS: lYork et al.ll20od IStoughton et"an 
[20021) . which provides well-calibrated, uniform optical images 
with a photometric accuracy of 2% -3%. The SDSS images 
are available in the ugriz filters, but the u and z images have 
very low S/N, and we concentrate our attention on the g, r, and 
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/ bands. We analyze the SDSS data following exactly the same 
procedures applied to the CGS. After registering the g and r 
images to the ; image, we extract isophotal intensity profiles as 
described in Section 4.1. As the SDSS images were observed 
in a drift-scan mode, they have both a very large field-of-view 
and a n exceptionally uniform background jPohlen & Truiilld 
120061) . This allo ws us to determine accurate sky val ues us- 
ing the method of iNoordermeer & van der Hulstl(l2007h . as de- 
scribed in Section 3. To convert the SDSS gri photometry into 
our stan dard BVRI system, we use the transformation equations 
given in lJester et aP (l2005h : 

B = g + 0.39(g-r) + 0.2l (9) 

V = g-0.59{g-r)-0.0l (10) 

7? = y-1.09(r- 0-0.22 (11) 

/ = /?-1.00(r- 0-0.21. (12) 

Figure[18]compares the surface brightness profiles from CGS 
with those derived from SDSS, for a subset of four relatively 
small galaxies (NGC 936, 1084, 1087 and 1090) that have well- 
measured sky values. The profiles are truncated at the radius 
where the isophotal intensity is 1 cr above the local sky value. It 
is clear that most of the profiles agree well with each other. The 
absolute value of the average profile differences is 0.24, 0.14, 
0.08, and 0.14 mag arcsec"- for B,V,R, and / bands, respec- 
tively. This level of discrepancy is not unexpected, given that 
all four of these galaxies were observed under non-photometric 
conditions in CGS, not to mention of additional uncertainties 
introduced by the photometric transformation from the SDSS 
to the CGS system. This comparison confirms that the basic 
reduction and calibration of the CGS data are sound. 

Despite the short exposures of the SDSS images (54 s), they 
have superior background uniformity and better sky determi- 



nation than CGS. These advantages translate to better sensitiv- 
ity in terms of surface brightness, by ~ 0.4, 0.2, 0.6, and 0.9 
mag arcsec"- in the B, V, R, and / bands, respectively. How- 
ever, the longer integration times, better seeing, and finer pixel 
scale of the CGS images imply that they have much higher S/N 
and sensitivity to compact structures compared to SDSS, typi- 
cally by a factor of ~ 4-5. 

11. SUMMARY 

We present a comprehensive isophotal analysis of optical 
(BVRI) images for a statistically complete, magnitude-limited 
sample of 605 bright, southern galaxies, observed as part of the 
Carnegie-Irvine Galaxy Survey (CGS). We discuss our strat- 
egy for determining the sky level and its error. It is challeng- 
ing to achieve very accurate sky subtraction with our images 
because the CGS galaxies are relatively large and the back- 
ground suffers from low-level non-uniformities due to residual 
flat-fielding errors. Nevertheless, crosschecks with internal and 
external data indicate that our calibration and sky subtraction 
strategies are robust, and that our quoted measurement uncer- 
tainties are sound. 

This paper focuses on the derivation of radial profiles of sur- 
face brightness, color, and various geometric parameters that 
characterize the shape and orientation of the isophotes. We 
construct composite brightness profiles as a function of Hub- 
ble type to highlight statistical trends. Non-exponential disks 
are seen in many SO and spiral galaxies. We perform a Fourier 
analysis of the isophotes to characterize their non-axisymmetric 
deviations from pure ellipses. The relative amplitude of the 
m = 1 mode effectively identifies lopsided structures in the light 
distribution, which we find to be common in our sample, espe- 
cially among late-type galaxies. Bars and spiral arms, by con- 
trast, are best revealed by the relative amplitude of the m = 2 
Fourier mode. We present a uniform set of quantitative mea- 
surements of bar size and bar strength, spiral arm strength, and 



TABLE 9: Parameter Differences for the Internal Comparison Sample 





AL125 


AD26.5 


Ae 


APA 


Am25 


Amtot 


Ai?20 


Aflso 


ARso 




(") 


(") 




n 


(mag) 


(mag) 


(") 


(") 


(") 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


Mean 


-1.21 


-3.86 


0.004 


5.53 


0.01 


0.01 


0.09 


0.25 


0.06 


Standard deviation 


20.42 


46.52 


0.059 


23.39 


0.09 


0.10 


0.81 


3.13 


11.83 



NOTE. — Mean difference and its standard deviation — Column 1: the isophotal diameter at 25 mag arcsec^^; Column 2: the 
isophotal diameter at 26.5 mag arcsec^ ; Column 3: ellipticity; Column 4: position angle; Column 5: total apparent magnitude 
within = 25 mag arcsec^^; Column 6: total apparent magnitude within the last reliable isophote (1 a above the sky); Column 
7: radius enclosing 20% of the total flux; Column 8: radius enclosing 50% of the total flux; Column 9: radius enclosing 80% of the 
total flux. 
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Fig. 15. — Internal comparison of the surface brightness profiles for the objects for which we have repeat observations. The data points are calculated by 
subtracting one profile from another. The horizontal dashed line in each panel is the weighted average of the residual data, and the two horizontal dotted lines mark 
the 1 (T photometric uncertainty of the two observations. 



lopsidedness amplitudes. 

Forthcoming papers will utilize the databases assembled here 
and in Paper 1 to explore a number of scientific issues, including 
the following. 

1 . Statistics of bars, bar properties, and their possible con- 
nection to spiral arms. 

2. Incidence of lopsidedness and its dependence on global 
galaxy and environmental parameters. 

3. Disk profiles, truncations, and correlations with color 
gradients. 

We thank the referee for a prompt and helpful review of this 



manuscript. This work was supported by the Carnegie Insti- 
tution for Science (L.C.H.), the UC Irvine School of Physi- 
cal Sciences (A.J.B.), the China Scholarship Council (Z.-Y.L.), 
and the Plaskett Fellowship of the Herzberg Institute of Astro- 
physics, National Research Council of Canada (C.Y.P). Z.-Y.L. 
is grateful to Professor X.-B. Wu of the Department of As- 
tronomy in Peking University for his support and helpful sug- 
gestions on this project. We made use of HyperLeda and the 
NASA/IPAC Extragalactic Database (NED), which is operated 
by the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics and Space 
Administration. Funding for the SDSS and SDSS-II has been 
provided by the Alfred P. Sloan Foundation, the Participating 
Institutions, the National Science Foundation, the U.S. Depart- 



24 



The Carnegie-Irvine Galaxy Survey. II. 



I 



0.5 


-0.5 
0.5 



o 

S -0.5 

CO 
(0 




-0.5 - 



0.5 - 



-0.5 - 



' 1 

B NGC 1688 



R 





T ' 1 ' r 

B NGC 1700 



R 





' 1 

B NGC 2196 




50 



100 



50 100 

r (arcsec) 




50 100 150 



Fig. 15. — continued 



ment of Energy, the National Aeronautics and Space Admin- 
istration, the Japanese Monbukagakusho, the Max Planck So- 
ciety, and the Higher Education Funding Council for England. 
The SDSS Web site is http : / /www . sdss . org[ 



Li et al. 



25 



0.5 


-0.5 

0.5 


I 

o 
<u 

-0.5 

O 
t-, 

CS 

tafl 

J. 

^ 

-0.5 

0.5 


-0.5 



1 ' 1 ' 1 ' 1 ' 
" B NGC 2310 '■. 


1 ' 1 ' 1 ' 
" B NGC 2434 


1 ' 


1 ' 1 ' 1 
" B NGC 2698 f 










T T 

-1 M 














III 

V 


1 

• > ! 


-| 1 1 1 1 1 — 

: ..A. 




H ' 1 ^ 

V 


T 








• i i 


n r( r < F 


1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 


H 1 1 ^ 


-- 


III 

" R 


1 


III 

" R 


1 


' R 










nr,r- = 








( t 


I ' ' ' 

I 




-| 1 1 1 1 1 1 h 

/ 


H ' 1 ^ 

I 










> 

<) 

(> 

1., , 






1,1,1,1, 


1,1,1, 


1,1, 


1 



50 100 150 50 100 150 

r (arcsec) 



50 



100 



Fig. 15. — continued 



26 



The Carnegie-Irvine Galaxy Survey. II. 




Fig. 15. — continued 



Li et al. 



CM 

I 



0.5 


-0.5 

0.5 




o 

S -0.5 

CO 



0.5 - 



-0.5 

0.5 


-0.5 



1 ' 1 ' 1 ' 
B NGC 7329 


- 




J»_ - 






-i — 1 — 1 — 1 — 1 — 1 — 

V 


h 

- 










-l 1 1 1 1 1 

R 


1 






!>- ^ 








-i — 1 — 1 — 1 — 1 — 1 


i> 

O- - 


h 

J> 




1,1,1 








1 ' 1 

B NGC 7590 



50 100 150 






T ' 1 ' r 

V NGC 1374 



H H 



H H 



B NGC 3672 



#-■' 



H H 



:*r*-i-.' 



J 1 L 



50 100 

r (arcsec) 



50 100 150 




Fig. 15. — continued 



28 



The Carnegie-Irvine Galaxy Survey. II. 




Fig. 16. — Internal comparison of derived parameters for objects with repeat observations. Observation 1 denotes values we adopt in the survey, while Observation 
2 gives the reference values for internal comparison. The dashed line denotes y = x. 
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HyperLeda 

Fig. 17. — External comparison of derived parameters between CGS and HyperLeda. The full sample available for comparison is shown as black squares, while 
the red open points mark the subset with the smallest measurement uncertainty. For Itot. Btot, B-V, and D25, the most reliable points are those that were observed 
under photometric conditions, that have well-measured sky values, and that do not suffer from contamination by nearby bright field stars. For e the red points 
highlight objects with minimal blight star contamination, and for the PA we additionally require that e > 0.3. The dashed fine denotes y = x. (A color version of this 
figure is available in the online journal.) 
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Fig. 18. — Comparison between the surface brightness profiles derived from CGS and SDSS. The SDSS data were analyzed in the same way as the CGS data, 
and their photometric system was transformed to ours as described in Section 10.3. Within each panel, the upper subpanel shows the two profiles, truncated at the 
radius where the isophotal intensity is 1 cr above the local sky background; the bottom subpanel shows the difference between the CGS and SDSS profiles. (A color 
version of this figure is available in the online journal.) 
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APPENDIX 

DATABASE OF ISOPHOTAL PARAMETERS 

Figures 19. 1-19.616 present the database of isophotal parameters for the 605 galaxies in CGS; we also include the 1 1 extra galaxies 
that are not part of the formal sample. One full-page figure is devoted to each galaxy, ordered sequentially following the numerical 
indices listed in Column 1 of Table 2. In the lower-left panel, we give the B, V, R, and / surface brightness profiles, followed by the 
B-I, V -I, and 7?-/ color profiles. The lower-right panel shows the radial profiles of e, PA, A3 (blue triangles) and B3 (red circles), 
A4 (blue triangles) and B4 (red circles), /1//0 (black triangles) and (f)\ (red circles), and h/h (black triangles) and 02 (red circles). 
The horizontal dashed lines in the e and PA subpanels denote the characteristic values of the galaxy. The top two panels show the BVI 
color composite image, displayed using an arcsinh stretch, and the star-cleaned /-band structure map, displayed using a linear stretch. 
These are helpful to show side-by-side with the isophotal plots to aid in the identification and interpretation of specific features in the 
1-D radial profiles. We only display three sample pages for illustration (Figures 19(a)-19(c); the full set of figures is available in the 
electronic version of the paper, as well as on the project Web site |http : / / cgs . obs . carnegiescience . edu| 
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Fig. 19(A). — Brightness profiles and isophotal parameters for the barred SBO/a galaxy NGC 1452. The color composite image is shown in the upper left, 
displayed using an arcsinh stretch, and the star-cleaned /-band structure map is shown in the upper right, displayed using a linear stretch. The lower-left panel gives 
the B, V, R, and / surface brightness profiles, followed by the B-I,V —I, and R— / color profiles. The lower-right panel shows the radial profiles of e, PA, A3 (blue 
triangles) and B3 (red circles), A4 (blue triangles) and B4 (red circles), Ii /Iq (black triangles) and 4'[ (red circles), and /2//0 (black triangles) and (^2 (red circles). 
The horizontal dashed lines in the e and PA subpanels denote the characteristic values of the galaxy. (A color version and the complete figure set (616 images) are 
available in the online journal.) 
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Fig. 19(B). — Brightness profiles and isophotal pai'ameters for the lopsided SBc galaxy NGC 3513. The color composite image is shown in the upper left, 
displayed using an ai'csinh stretch, and the stai'-cleaned /-band structure map is shown in the upper right, displayed using a linear stretch. The lower-left panel gives 
the B, V. R, and / surface brightness profiles, followed by the B-I,V —I, and R— / color profiles. The lower-right panel shows the radial profiles of e, PA, A3 (blue 
triangles) and B3 (red circles), A4 (blue triangles) and B4 (red circles), /| /Iq (black triangles) and 4'[ (red circles), and /2//0 (black triangles) and (^2 (red circles). 
The horizontal dashed lines in the e and PA subpanels denote the characteristic values of the galaxy. (A color version of this figure is available in the online journal.) 
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Fig. 19(C). — Brightness profiles and isophotal parameters for the ban'ed SBb galaxy NGC 7329. The color composite image is shown in the upper left, displayed 
using an arcsinh stretch, and the star-cleaned /-band structure map is shown in the upper right, displayed using a lineal' stretch. The lower-left panel gives the B, V, R, 
and / surface brightness profiles, followed by the B—I, V—I, andR—I color profiles. The lower-right panel shows the radial profiles of e, PA, A3 (blue triangles) and 
B3 (red circles), A4 (blue triangles) and B4 (red circles), I[ /Iq (black triangles) and (red circles), and L/Iq (black triangles) and (/>2 (red circles). The horizontal 
dashed lines in the e and PA subpanels denote the characteristic values of the galaxy. (A color version of this figure is available in the online journal.) 



